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A central issue in the study of polymer physics is to understand the relation between the geometrical properties 
of macromolecules and various dynamics, most of which are encoded in the Laplacian spectra of a related graph 
describing the macrostructural structure. In this paper, we introduce a family of treelike polymer networks 
with a parameter, which has the same size as the Vicsek fractals modeling regular hyperbranched polymers. 
We study some relevant properties of the networks and show that they have an exponentially decaying 
degree distribution and exhibit the small-world behavior. We then study the Laplacian eigenvalues and 
their corresponding eigenvectors of the networks under consideration, with both quantities being determined 
through the recursive relations deduced from the network structure. Using the obtained recursive relations 
we can find all the eigenvalues and eigenvectors for the networks with any size. Finally, as some applications, 
we use the eigenvalues to study analytically or semi-analytically three dynamical processes occurring in the 
networks, including random walks, relaxation dynamics in the framework of generalized Gaussian structure, 
as well as the fluorescence depolarization under quasiresonant energy transfer. Moreover, we compare the 
results with those corresponding to Vicsek fractals, and show that the dynamics differ greatly for the two 
network families, which thus enables us to distinguish between them. 



PACS numbers: 36.20.-r, 64.60. aq, 89.75.Fb, 05.40.Fb 



I. INTRODUCTION 

A fundamental issue in the study of complex systems is 
to unveil how the structural properties affect various dy- 
namics, many of which are related to the exact knowledge 
of the eigenvalues and eigenvectors of Laplacian matrix. 
Examples include relaxation dynamic in the framework 
of generalized Gaussian structure (GGS) 1 -, fluorescence 
depolarization by quasiresonant energy transfer^* 3 -, stan- 
dard discrete-time random walks^, and continuous-time 
quantum walks^S, and so on. In addition to dynami- 
cal processes, Laplacian eigenvalues and eigenvectors are 
also relevant to diverse structural aspects of complex sys- 
tems, such as spanning freest and resistance distance 8 -. 
Thus, it of theoretical interest and practical importance 
to derive exact analytical expressions of Laplacian eigen- 
values and eigenvectors for complex systems, which can 
lead to extensive insights in the contexts of topologies 
and dynamics. 

Given the wide range of applicability, the study of 
Laplacian eigenvalues and eigenvectors has been sub- 
ject of considerable research endeavor for the past few 
decades. Thus far, the Laplacian eigenvalues for some 
classes of graphs have been determined exactly, includ- 
ing regular hypercubic lattices^, dual Sierpinski gas- 
ket o 10 ' 11 , Vicsek fractal s 12 ' 13 , dendrimer also known as 
Cayley treeM, and Husimi cact i 15 ' 16 . Recent empirical re- 
search indicated that some real- life networks (e.g., power 
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grid) display small-world behavio r 17 ' 18 . Moreover, these 
networks are simultaneously characterized by an expo- 
nentially decaying degree distribution^ 8 , which cannot 
be described by above-mentioned networks. However, 
related work about Laplacian eigenvalues and eigenvec- 
tors for small-world exponential networks is much less, 
notwithstanding the ubiquitous nature of such systems. 

In this paper, we define a category of treelike polymer 
networks controlled by a parameter, which is built in an 
iterative way. The networks have the same size as that of 
Vicsek fractal a 19 ' 20 corresponding to the same parameter 
and iteration. According to the construction, we study 
some structural properties of the networks, showing that 
they have an exponentially decaying degree distribution, 
and display the small-world property. Moreover, the net- 
works can be assortative, uncorrelated, or disassortative, 
relying on the parameter. Then, by applying the tech- 
nique of graph theory and an algebraic iterative proce- 
dure, we study the Laplacian eigenvalues and eigenvec- 
tors of the networks, obtaining recursive relations for the 
eigenvalues and eigenvectors, which allow for determining 
exactly the full eigenvalues and eigenvectors of networks 
of arbitrary size. 

In the second part of this work, by making use of the 
obtained Laplacian eigenvalues, we study three classic 
dynamics for the small-world polymer networks, such as 
trapping with a single trap, relaxation dynamics in the 
GGS framework, and the fluorescence depolarization un- 
der quasiresonant energy transfer. For the trapping prob- 
lem, we study two particular cases: in the first case the 
trap is fixed at the central node, while in the other case 
the trap is distributed uniformly. For both cases, we de- 
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rive explicit formulas for the average trapping time and 
obtain their leading scalings, which follow different be- 
haviors, showing that the position of trap has a substan- 
tial effect on the trapping efficiency. For the GGS, we 
determine three interesting quantities related to the re- 
laxation dynamics, i.e., the averaged monomer displace- 
ment, storage module and loss module. Finally, we dis- 
play the behavior of the fluorescence depolarization. For 
the three dynamics, we also present a comparison for the 
behaviors between the small- world polymer networks and 
Vicsek fractals, and show that they differ strongly. 



II. NETWORK CONSTRUCTION AND PROPERTIES 

In this section, we first introduce a family of treelike 
small-world polymer networks with an exponential de- 
gree distribution, then we study some relevant properties 
of the networks. 



A. Construction method 

The networks being studied have a treelike structure, 
and are constructed in a deterministically iterative way. 
Let U g (g > 0) denote the networks after g iterations. 
For g — 0, Uo consists of an isolated node, called the 
central node. For g = 1, / (/ is a positive integer) new 
nodes are generated connecting the central node to form 
Ui. For g > 1, U g is obtained from U g ~\ by attaching 
/ new nodes to each node in U g -\. Figure [TJ illustrates 
schematically the first several iterative construction pro- 
cesses of a particular network for the case of / = 3. 



FIG. 1. (Color online) Construction of a special network cor- 
responding to / = 3. 

According to the construction approach, it is easy to 
derive that at each iterative step gi (gi > I), the number 
of newly generated nodes is L(gi) = /(/ + l) 9i_1 . Then 



the total number of nodes at each generation g is 

N g = i + J2L( 9i ) = (f + iy, 



(1) 



and the total number of edges in U g is E g = N g — 1 = 

if + iy-i. 

In fact, the networks being studied are self-similar, 
which can be seen from another construction approach. 
As will be shown below, the central node of U g has the 
largest degree, we thus also call it hub node. Let h g de- 
note the central node of U g . Then, U g can be constructed 
alternatively as follows, highlighting its self-similarity, see 
Fig. [2j To generate U g , we create / + I replicas of U g -i, 

and label them as U { g %, uf\, U { g 2 \,.. 

tively. Moreover, let h^ c } 1 (x — 0, 1, 2, . 



( f) 

, Uy_\, respec- 
. , /) denote the 



hub of the U ( g x \. Then, for each (i — 1,2,..., /), 

we introduce an additional edge connecting its hub node 

Thus, through the two steps of 



(i) 



to the node h g _ 



replication and connection, we obtain U g with h g _-y being 
its hub. 




FIG. 2. (Color online) Second construction method of the 
small- world polymer networks. U g consists of / + 1 copies 



of U, 



AO) 



9 -i, denoted by U*_ u ^ g _ t , ^ s _ lr . .,^ 9 _i, 



U, 



(i) 



v. 



(2) 



■ ■ ;U, 



(/) 



which are 



connected to each other to form U g 
between the central node of each I/g!^ 



by adding a new edge 
(* = 1, 2, . . . , /) and 



the central node of U 



(0) 



Note that the numbers of nodes and edges of the net- 
works under consideration are identical to those cor- 
responding to Vicsek fractal s 19 ' 20 , but their structural 
properties differ greatly from those of Vicsek fractals, as 
we will show. 



B. Structural properties 



We proceed to present some important structural prop- 
erties of U gi including degree distribution, average path 
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length, diameter, and degree correlations. 

1. Degree distribution 

For a network, its degree distribution P(k) is defined as 
the probability that a randomly chosen node has a degree 
of k. Let ki(g) be the degree of node i in U g . Assume that 
node i entered the networks at generation gi (<?j > 0), 
then ki(gt) = 1. By construction, at each subsequent 
iteration, / new nodes will be generated linking to node 
i. Thus, the degree of node i evolves as 

ki(g) = ki(g-l) + f. (2) 

Considering ki{gi) — 1, Eq. ([2]) is solved to yield 

fciG?) = ! + /(<?-<?*), (3) 

which provides the degrees of all nodes except the central 
one. We label the initial central node by 0; then the 
degree of node in U g is 



ko{g) = fg, 



(4) 



which is the highest among all nodes. 

Equations §5§ and dU show that the degree spectrum 
of U g is discrete and that all nodes generated at the same 
generation have the same degree. Thus, in U g , the num- 
ber of possible node degrees is g+1, which is in sharp con- 
trast to that for Vicske fractals, where only three types of 
degrees exist, that is, 1, 2 and /. It follows that the cu- 
mulative degree distribution 21 of the networks addressed 
is given by 



oo 

E 

k'=k 



Using Eq. ©, 
P(g'<<t> = g- 



we have P r , 
fc-i 
/ 



P(k') 



(•5) 



Pr, 



>(*) = E 



. Hence, 
W) (/ + !)*" 



N„ 



(f+iy 



k—l 

f 



(6) 

which decays exponentially with k. It is the same with 
degree distribution P(k), see^i for explanation. 



2. Average path length 

The average path length represents the average of 
length of the shortest path between two nodes over all 
node pairs. Assume that each edge in U g has a unit 
length. Then the length of the shortest path between 
nodes i and j in U g , denoted by dij(g), is the minimum 
length for the path connecting the two nodes. Let d g 
represent the average path length of U g , defined by: 



Stot(g) 



N g (N g - l)/2 



(7) 



where Stot(s) is the sum of dij(g) over all pairs of nodes, 
i.e., 



Statig) =J2d ij (g). 



(8) 



We note that in Eq. ([8]) , for a pair of nodes i and j (i ^ j ) , 
we only count dij(g) or dji(g), not both. 

Let & g and <d g the sets of nodes generated at iteration 
g or earlier, respectively. Then Stot{g) can be recast as 

Stot{g)= Yl d v(9)+ Yl d v(9)+ Y d ^(9), 
ie& g ,jee g iee g ,jee g iee g ,jee g 

(9) 

It is evident that the third term on the right-hand side 
(rhs) of Eq. ^ is exactly 5 to t(.9 — l) f i-e., 



Y2 d ij(d) = S to t(g- 1) 
iee g , jee g 



(10) 



For the first two terms on the rhs of Eq. ©, according 
to the first network construction method, they can be 
evaluated as 

J2 d l3 {g) = / [{N g ^f + 2S tot (g - 1)] , (11) 

i£© s , j£@ g 

and 

E d *M = / 2 Stot(5 - 1) + / AW/ N g -! - 1) , 
»ee 9 , jee g 

(12) 

respectively. 

Plugging Eqs. (jlOH12|) into Eq. © leads to 

Stot(g) = (/+ l) 2 S tot (g - 1) + / (/ + l)(iV 9 ~i) 2 - /JV S _! 

9-1 

= (/ + i) 2s Stot(o) + / (/ + [(/ + 1) 2 Cs - 1 " i) (iv i ) 2 ] 

i=0 

-/E[(/ + 1 ) 2(9 " 1 " I)jV J ( 13 ) 



Substituting 5 tot (0) = and N t = (/ + 1) 4 into Eq. pB . 
we can obtain the exact expression for Stot(g) as 

Stat GO = if 9 - + If 9 - 1 + (/ + I) 5 " 1 • (14) 
Inserting Eq. (fbl)) into Eq. (JTJ) gives 

j (fg-W + lf^ + if + l) 3 - 1 



(/ + l) 9 [(/ + l) 9 -l]/2 
_ 2(/g-l)(/ + l) 9 + 2 

(/ + l) 9+1 -(/ + l) ' 



(15) 



Recalling N g = (/ + l) s as given in Eq. ([T]), we have 
g = In 2V fl / ln(/ + 1), both of which enable us to write d g 



in term of network size N g as 

- _ 2(/ lnA g /ln(/ + l)-l)A g + 2 
(f + l)N g -(/ + 1) 
2/ A s lnA g 2 1 



(/ + l)ln(/ + l) A g -1 / + 1A S -1 



.(16) 
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When the network size is large enough, we have 
t ~ 2/ 



(1 + /) ln(l + /) 



lnN g , 



(17) 



which increases logarithmically with the network size g, 
showing that the networks display the small-world be- 
havior—. 



3. Diameter 

We have shown that the treelike polymer networks are 
small-world, since their average path length grows as a 
logarithmic function of network size. In addition to av- 
erage path length, sometimes, diameter is also used to 
characterize the small- world phenomenon, since small di- 
ameter is consistent with the concept of small- world. For 
a network, its diameter is defined as the maximum of the 
shortest distances between all pairs of nodes in the net- 
work. Let diam([/g) denote the diameter of U g , below we 
will compute analytically diam(t/ ff ) and show that it also 
scales logarithmically with the network size. 

Clearly, at step g = 1, diam(L/i) equals 2. At each 
iteration g > 1, we call newly-generated nodes at this it- 
eration active nodes. Since all active nodes are connected 
to those nodes existing in U g —i, it is easy to see that the 
maximum distance between an arbitrary active node and 
those nodes in f7 ff _i is not more than diam({/i) + 1 and 
that the maximum distance between any pair of active 
nodes is at most diam(£/i) + 2. Hence, at any iteration, 
the diameter of the network increases by 2 at most. Then 
we get 2(g + 1) as the diameter of U g , which is equal to 
2(logy +1 N g + 1) growing logarithmically with the net- 
work size. This again indicates that the networks under 
study are small-world. 



is in the interval [—1,1]. If the network is uncorrelated, 
r equals zero. Disassortative networks have r < 0, while 
assortative graphs have r > 0. Let r(f,g) be the Pearson 
degree correlation coefficient of U g . By definition, r{f,g) 
is given by 



r(f,g) 



1 2 



. (18) 



where where ji and fc, are the degrees of the nodes at the 
two ends of the ith edge in U g , where i 6 {1, 2, . . . , E g }. 

The three terms in numerator and denominator in 
Eq. (fT8|) can be evaluated as 



J2 3i ki = (3/ + 7) (/ + 1 Y - 2/ V - 7fg - 3/ - 7 , (19) 
£ = -^(/+5)(/ + l) 9 + 2/<?+^(/ + 5) , (20) 



and 



-i(3/ 2 + 15/)-i(/ 2 + 9/ + 16), 

(21) 



4. Degree correlations 

For a network, its degree correlations^ 2 , can be de- 
scribed by the Pearson correlation coefficient r, which 



respectively. Inserting Eqs. ([T9" l) -(|2l] ) and E g = (/ + l) s - 
we can arrive at the explicit expression 



1 into Eq. ([T 
for r(f,g) as 



r(f,9) 



V[(f + I) 3 - IE/ + l) g + 6g - 1] - / 2 [(/ + iy - l] 2 + 3[(/ + 1)9 - l] 2 + 8g[(f + 1)9 (g - 1) + g + 1] 
(/ + 1) 29+1 (/ + 7) - 2/(/ + l)B[f + 8 - (/ + 5)g + 6/.g 2 ] - 14(/ + 1)9 + /[/ + 8 - 2(/ + 5)g - W] + 7 ' 

(22) 

I 



In Fig. [3l we report the exact result for r(/, g) provided 
by Eq. (g21). From Fig. H it is obvious that for / = 1, 2, 
r(f, g) is positive; for / = 3, r(/, g) equals zero; while for 
/ > 4, r(f,g) is negative. 



Equation (|22|l shows that for very large g, we have 
(/ + l) 29 (-/ 2 + 2/ + 3) 



r{f,g) 



(/ + l) 2s (/ 2 + 8/ + 7) 
/-3 

f + r 



(23) 



which decreases with /. When / = 1 and / = 2, r(f,g) 
is equal to j and ^, respectively. Thus, for / = 1 and 
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FIG. 3. (Color online) Pearson correlation coefficient r(f,g) 
of U g as a function of / and g. 



A. Eigenvalues 



Let A g = [-A»j](/+i)ffx(/+i)« denote the adjacency ma- 
trix of Ug , where Aij — Aji — 1 if nodes i and j are ad- 
jacent, Aij = Aji — otherwise, then the degree of node 
i is di = J2jeu g Lct D s = diag(di, d 2 , 
denote the diagonal degree matrix of J7 ff , then the Lapla- 
cian matrix of U g is defined by L s = D g — A g . 

We first study the eigenvalues of U g , leaving the eigen- 
vectors to Subsection IIIIBI By construction, it is easy to 
see that A g and D g obey the following relations: 



A„ = 













h-1 \ 





/ 



(24) 



/ = 2, U g is assortative. When / = 3, r(f,g) is equal to 
0, indicating that the network is uncorrelated. When / 6 
[4, oo), r(f,g) is negative. Concretely, when / increases 
from 4 to oo, r(f,g) decreases from - jj to —1, showing 



that U g is disassortative. 



and 



D, 








I 






9-1 



V 









Ia-1 / 



(25) 



The phenomenon that the Pearson degree correlation 
coefficient r(f,g) decreases with / can be explained 
heuristically as follows. Note that there are E g = 

(/ + l) 9 - 1 edges in U g , which means that for those in which each block is a (/ + l) 9 " 1 x (/ + I) 9 " 1 matrix 



Ng-i = (/ + 1 ) 9 1 old nodes having a degree higher than 
one, they have 2E g - X +L(g) = (/+2)(/ + l)f- 1 — 2 neigh- 
boring nodes, among which L{g) = /(/ + 1) 9_1 neighbors 
are those newly generated nodes with a single degree. 
Thus, for large g, the fraction of neighbors with single 
degree is approximatively equal to //( / + 2), which is an 
increasing function of /, meaning that in networks cor- 
responding to larger /, the average degree of neighbors 
of old nodes is smaller. 



III. LAPLACIAN EIGENVALUES AND THEIR 
CORRESPONDING EIGENVECTORS 

Although for general graphs, it is a challenge to de- 
termine their Laplacian eigenvalues and eigenvectors, as 
will be shown, for U g this problem can be settled. 



and I g -i is the (/ + l) 9 " 1 x (/ + l) 9 " 1 identity matrix. 
Thus, the Laplacian matrix of U g satisfies the following 
recursive relation: 



/ L g _l + flg-l —Ig-1 —Ig-1 



-I 9 -l 



-I 9 -l 



I 9 -l 



1,-1 



"I 9 -l \ 






(26) 



J 



Obviously, the problem of determining Laplacian 
eigenvalues of U g is equivalent to finding the roots of 
characteristic polynomial P g (X) of L s . To find the eigen- 
values of L g , we just need to determine the roots of P g ( A) , 
which reads: 
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P t (A)=det(AI ff -L g ) 

/ (A - /)!,_! - L s _! Vi 



del l 9~l 



= {dct[(A- - dot 



(A - l)Vi 




= {det[(A - l)I g -i]} f ■ det 



o 

(A - l)Vi 



1 T T 

A-l 9-1 L 9- 

1 T 
1 T 





I9-1 \, 


















• (A 


- 1) Vi J 




9-1 


••• I9-1 \ 




U 


■ • • u 




9-1 


n 







■•• 19-J 







••• 


\ 


I9-1 


■•■ 








Vi ••• 






1 T 

— A 9-l 







(27) 



where we have used the elementary operations of matrix. 
Based on the results ir>22, -Pg(A) can be expressed as 



P g (A) = {det[(A-l)I g _ 1 ]} / det 



X — f — 



A- 1 



Ig-l - L g _l 

(28) 



Hence, -Pg(A) can be further recast recursively as 
F 9 (A) = (A-l) /(/+1)S > 9 -i(^(A)) 



(29) 



where ip(X) = A — / — t^t- This recursion relation 
provided in Eq. (|29]l is very useful for determining the 
eigenvalues and eigenvectors of the Laplacian matrix for 
U g . Note that P 9 _i(A) is a monic polynomial of degree 
(/ + then the exponent of j^j in P g _i(ip(X)) is 

(/ + 1) 9_1 , and the exponent of factor (A — 1) in P g (X) is 



/(/ + l)^ 1 - (/ + l)^ 1 = (/ - !)(/ + 1)^-1 _ (30) 

Therefore, U g has Laplacian eigenvalue 1 with multiplic- 

ity (/-l)(/ + l) fl - x . 

It is evident that U g has (/ + I) 9 Laplacian eigen- 
values, denoted by Af , Af, • ■ • , X^ +1 -. g , the set of which 

is represented by A g , i.e., A g — {Af, Af, . . . , A^ +1 ^ s }. 
In addition, without loss of generality, we assume that 

On the basis of above analysis, 

such 

Ag 1 '' contains all eigenvalues equal 
to 1, while A'if' includes the remain eigenvalues. Thus, 

A« = {1,1,1,..., 1,1}, (31) 
" , ' 

(/-1X/+1) 9 - 1 

where the distinctness of elements is neglected. 



Af < Af < . 



< \ 9 
- A (/+i) 9 



A g can be divided into two subsets Ag 1 ^ and A^ 
as Ag = A 9 1} U A 9 2) . 

,(2) 



The remaining 2(f + l) 9 1 eigenvalues belonging to 

(2) 

Ag are determined by P g _i(<^(A)) = 0. Let the 
2(f + l) 9 ^ 1 eigenvalues be Af , Af, ■ • ■ , Afy +1 N/_i, respec- 



tively. That is, Ag 2) = {Af , Af , 



, Xi, . ,_ 1)a _i}. For con- 



venience, we assume that A? < Af < ... < X 9 ±. 
Equation (|2T)|) shows that for any element in A ff _i, say 
X 9 ~ G Ag_i, both solutions of A — / — -J^ = Af -1 are 



(2) 

in Ag ' . It is clear that A — / • 
to 



/ 

A-l 



A 



9-1 



is equivalent 



A 2 - (A?" 1 + / + 1)A + A?" 1 = 0, 



(32) 



the two roots of which are denoted, respectively, by Af 
and -Af+(/+i).t7-i ' since these notations give a natural in- 
creasing order of the eigenvalues of U g , as will be shown 
below. 

Solving the quadratic equation provided by Eq. (|3"2"1) . 

we obtain the two roots to be Af = ri(Af _1 ) and 

Af +(/+ i )9 -i = 7- 2 (Af _1 ), where r^Af" 1 ) and r 2 (Af _1 ) 
are 



^(Af- 1 ) = \ ( Af" 1 + / + 1 - ^(Af- 1 + / + l) 2 - 4AP 1 



and 



(33) 



r 2 (xr 1 ) = \ (V 1 + / + 1 + \/(Ar 1 + /+i) 2 -4Ar 1 

J34) 

respectively. Thus, in this way each eigenvalue Af in 

(2*) 

Ag_i gives rise to two new eigenvalues in A g . Inserting 
each Laplacian eigenvalue of U g -\ into Eqs. (|33|) and ([34]) 

(2) 

generates all the elements of Ag . Considering the ini- 
tial value A = {0}, by recursively applying Eqs. (|3"3")l 
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and ([34]) . the Laplacian eigenvalues of U g can be fully 
determined. 

It is easy to prove that the two roots, r 1 (X 9 ^ 1 ) and 
?"2(Af~ ), of Eq. (|32j) monotonously increase with 
and both lie in intervals [0, 1) and (1, +oo), respectively. 
Thus, for any eigenvalue in Af~ 6 A g _i, r\{X 9 i ~ ) < 1 < 
?"2(Af _1 ) always holds. In addition, the following conclu- 
sion can be reached based on simple argument. Assum- 
ing that E g -x = {Af~^Af~V..,A?7* .J, then A g 2) 



can be generated via Eqs. 1|33|) and ([34 



{A? 



^2) • 



1 < A ?/ . 



\ 9 

'•' /V 2(/+l)9 



} satisfying Af < A| 



that is, A ff } 
< X 9 

- (/+i) 3 



< 



that Ag 1 ' contains (/ — !)(/ • 



-i - /K u+iy- 1 +2 



< X s , 



Recall 



2(/+l)"- 1 ' 

elements 1, we now 



have gotten the whole set of Laplacian eigenvalues for U g 
to be A fl =A S 1) UA S 2) . 

In order to see the distribution of the Laplacian eigen- 
values for U g . We use Eqs. (f3"3"]l and (IMl) to determine 
the eigenvalues of a specifical network corresponding to 
/ = 4 and g — 5. In addition, by diagonalizing the 
associated Laplacian matrix, we also compute numeri- 
cally the eigenvalues and their multiplicities, which are 
in complete agreement with those analytical results, con- 
firming that the theoretic approach is valid. In Fig.@|a), 
we display as a histogram, for the result of the network 
corresponding to / = 4 and g = 5, thus having a size 
N$ = 3125. Furthermore, we also present in Fig. Sth) 
the histogram for the corresponding Vicsek fractals with 
/ = 4 and g = 5. 

By comparing Figs. Ufa) and (b), we can see that num- 
ber of distinct eigenvalues in the small-world network is 
much less than its corresponding Vicsek fractal. Note 
that in U g , the distinct degree values for nodes are g + 1, 
while for corresponding Vicsek fractals, the degree values 
are 3 (all node have degree 1, 2, or /). The reasons for 
the interesting phenomenon that Vicsek fractals display 
a larger heterogeneity in the Laplacian spectrum but a 
far smaller heterogeneity in the degree values deserves 
further study in the future. In addition to the number 
of dissimilar eigenvalues, the difference of eigenvalues are 
also obvious for these two networks. For instance, the 
maximum eigenvalue, Ag lax , of the small- world polymer 
network is substantially higher than that of the Viscek 
fractal. As we will show, these differences of Laplacian 
spectra between the two networks will lead to different 
behaviors for various dynamics taking place on them. 



B. Eigenvectors 
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FIG. 4. Number of distinct eigenvalues for a small-world poly- 
mer network (a) and its corresponding Vicsek fracal, with 
/ = 4 and g = 5 for both networks. 



(\I g — L g )v = 0. We distinguish two cases: A G Ag 1 - 1 and 
A G A g , which will be separately treated as follows. 
For the case of A G Ag 1 ^, in which all A = 1, equation 



(AI S — Lg)v = becomes 



/(l-/)I 9 _i-L s _ : Iff - 1 

Ig-l 

Ig-1 



I 



9-1 








/ 



l>2 
V 3 



= 0, 



V / 

(35) 

where vector Vi (1 < i < f + 1) arc components of v. 
Equation (|35|) leads to the following equations: 



Analogous to the eigenvalues, the eigenvectors of L g 
can also be derived directly from those of L s _i. Assume 
that A is an eigenvalue of Laplacian matrix for U g , the 
corresponding eigenvector of which is v G +1 ) , where 
p,(/+!) 9 is the (/ + l) 9 -dimensional vector space. Then 
the eigenvector v can be determined by solving equation 



V2 + V 3 



vi = 0, 
Vf+i = 0. 



(36) 
(37) 



In Eq. (|36j) . v\ is a zero vector. Let Vi = 
(wiji, Uj,2j ■ ■ ■ ! v i,(f+i)s) j then, Eq. (l37j) is equivalent to 
the following equations: 



8 



V2,l + «3,1 + ... + W/+l,l = 

V2,2 + V 3>2 + ... + Vf +1) 2 = Q 



' ; 2,(/+l)s-! + v 3,(f+1-) ! '- 1 + ••• + «/+l,(/+l)s-l — 

The set of all solutions to any of the above equations 
consists of vectors of the following form 



v 4.j 
V Vf+ij J 





V o J 



«2,j 



1 



k-m—l, 



o 

V i J 



(38) 



where kij , k 2 j , ■■■ , ar e arbitrary real numbers. In Eq. (|3"8"j). the solutions for all the vectors Vi (2 < i < f + 1) 

can be rewritten as 



\ V l+1 I 



( 



-1 -1 

1 
1 


















&2,2 




&3.1 


&3,2 



V ••• 1 / 



■•• h,(f+i)9-i 



(39) 



where fc^- (l < i < / - 1; 1 < j < (/ + l) s_1 ) are 
arbitrary real numbers. Using Eq. (131)1) . we can obtain 
the eigenvector v associated with the eigenvalue 1. Fur- 
thermore, we can easily check that the dimension of the 



eigenspace of matrix L 9 corresponding to eigenvalue 1 is 
(/-!)(/ + I) 3 " 1 - 

We proceed to address the case of A G A^'. For this 
case, equation (AI g — L g )t; = can be rewritten as 



/ (A - f)l g _ x - L s _i 

Vl 



(A-1)I S 




i 

(A - 1)V 



l 9-l 






/ »1 \ 



V-2 
V 3 



(X-l)l a -J \v f+ J 



(40) 



where vector v% (1 < i < f + 1) are components of t;. 
Equation (|4"0)) leads to the following equations: 

[(A - f)I g _ t - L a _i] wi + v 2 + ■ ■ ■ + v f+1 = 0, (41) 
vi + (X — l)vi = (2 <« </ + !). (42) 



Resolving Eq. (|42|) yields 
1 



A- 1 



vi (2 <»</ + !) 



Inserting Eq. (1431) into Eq. (|4"Tj) results in 



A-/- 



/ 

A- 1 



lg-1 



■>g-\ 



(43) 



wi = 0, (44) 



which indicates that v\ is the solution of Eq. (|4"T|) while 
(2 < i < f + 1) are completely determined by Vi 
via Eq. (|4"3")l . As demonstrated in Eq. (|2^|) . if A is an 
eigenvalue of L g , then <ys(A) = A — /— -j-rr is an eigenvalue 
of lig-i- Thus, Eqs. (|4"4"|) and (|2§f implies that Ui is an 
eigenvector of L s _i corresponding to eigenvalue A — / — 
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jry, while 



Let F sum (g) denote the sum term on the rhs of Eq. (|46l) . 



v 3 



(45) 



is an eigenvector of L g associated with eigenvalue A. 

Since for the initial graph Uq, its Laplacian matrix Lo 
has only one eigenvalue with corresponding eigenvector 
(1) T ; by recursively applying the above process, we can 

(2) 

obtain all the eigenvectors corresponding to A 6 A g . 

In this way, we have completely determined all eigen- 
values and their corresponding eigenvectors of U g . In 
the following text, we will use these obtained results, es- 
pecially those for eigenvalues, to study some dynamical 
processes taking places in U g , including random walks 
with a trap, relaxation dynamics in the GGS framework, 
and depolarization of fluorescence by Foster quasireso- 
nant energy transfer. 



IV. TRAPPING PROCESS 

In this section, we study trapping problem in the small- 
world polymer networks. The trapping problem is a par- 
ticular kind of random walks with a trap fixed at a posi- 
tion, absorbing all particles visiting it. In the process of 
random walks, at each time step, the particle (walker), 
starting from its current location, moves to any of its 
nearest neighbors with equal probability. One of the pri- 
mary quantities related to trapping problem is trapping 
time (TT)21 The TT for a node is defined as the mean 
first-passage time (MFPT) for a particle starting from 
the node to the trap. Let Fij(g) denote the MFPT from 
node i to node j. Below we will focus on two cases of 
trapping problem. In the first case, the trap is fixed on 
the central node, while in the other case, the trap is uni- 
formly distributed over the whole networks. 



A. Trapping with a trap fixed on the central node 

We first consider the case of trapping in U g with the 
perfect trap being located at the central hub node h g . 
In this case, the quantity we are concerned with is the 
average trapping time (ATT), Fh(g), which is the average 
of F it h g (g) over all possible starting points in U g . That 



1 Ng 

F h (g) = w Y, F ^9)- 



(46) 



We next study analytically Fh(g) by using the second 
construction method of the networks, showing how Fh{g) 
changes with the network size N g . 



i.e. 



Then, 



ieu„ 



F h (g) = 



Fsum(g) 



(47) 



(48) 



Thus, we reduce the problem of determining Fh(g) to 
evaluating F sum (g). To find F snm (g), we should deter- 
mine some intermediary quantities. First, for all g > 0, 
Fh g ,h g {g) — 0. On the other hand, according to the pre- 
vious results obtained by various technique s] 25 ' 26 , we have 

F H?*.(S) = F h ^i9) = 2iV ff _ x - 1 = 2(/ + l)*- 1 - 1 

. (49) 

for all 1 < i < /. Then, from the second construction of 
the networks, we obtain 



F SU xn(g) 

f r 

= £ *i,*„(ff)+£ £ [ F i, h pi0) +F Hp,hpW 

= F sum {g - 1) + f[F sum {g - 1) + JV fl _i(2JV fl _i - 1)] 

= (/ + i)F sum (. 9 - 1) + /(/ + ly-'m + ly- 1 - 1] . 

(50) 

Considering F SU m(0) = 0, Eq. ([50]) is solved to yield 
F sum (g) - 2(/ + l) 2 ^ 1 - (/ + \y-\fg + 2) . (51) 



Substituting Eq. (|51|) into Eq. (|48|) . we arrive at the 
closed- form expression of Fh{g) as 



F h (g) =2(/ + l) 



9-1 



/g + 2 
/ + 1 



(52) 



We next show how to represent Fh{g) in terms of the 
network size N g , with a goal to obtain the relation be- 
tween these two quantities. Recalling Eq. ([T]), we have 
g = In N g / ln(/ + 1), which enables to write Fh(g) in the 
following form: 



F h (g) = 



2N n 



f + 1 (/ + 1) m(/ + 1) / + 1 



(53) 



Equation (|53j) provides an explicit dependence relation 
of Fh(g) on N g and parameter /. For a sufficiently large 
system, i.e., N g —> oo, the dominating term of Fh{g) is 



F h (g)^ 



/ + !' 



(54) 



which increases linearly with the system size. This linear 
scaling of ATT on the network size is in sharp contrast 
to the superlinear scaling of ATT in Vicsek fractals with 
the central node as the trap 27 ' 28 . 



10 



B. Trapping with the trap uniformly distributed 



In Subsection IIV A[ we have discussed the trapping 
problem in U g with an immobile trap positioned at the 
central node. Here we study another case of trapping 
problem in U g with the trap uniformly distributed over 
the whole networks. In this case, we are concerned with 
the quantity F g defined as the average of MFPT (g) 
over all pairs of source point i and target point j in the 
networks: 



F„ 



1 



1=1 j=l 



(55) 



Let F to t(g) denote the summation term on the rhs of 
Eq. ' 



Then, 



--1 3 = 1 



Ftot(g) 



(56) 



(57) 



which is actually the ATT when the trap is uniformly 
distributed. Notice that the quantity F g involves a dou- 
ble average: the first one is over all the source points to a 
given trap, the second one is the average of the first one. 

In order to compute F g , we use the relation governing 
resistance distance and MFPTs between two nodes in a 
connected grap h 29 ' 30 . For this purpose, we look on U g 
as an electrical network^ by considering each edge in 
U g to be a unit resistor—. Let Rij(g) be the effective 
resistance between two nodes i and j in the electrical 
network corresponding to U g . Then, the following exact 
relation 

F ij (g) + F ji (g)=2E g R ij (g) (58) 
hold o 29 ' 30 , and Eq. (|55| can be recast as 



F tot (g) = E g J2J2 Ri ^)- 

i=l j=l 



(59) 



Applying the previous result o 33 i 34 , the sum term of ef- 
fective resistance between all pairs of nodes in U g can be 
evaluated as 



(60) 



i=l j=l 

Then, Eq. (|5"5)) becomes 



Ng 1 



(61) 



Having expressing F g in terms of the sum of the recip- 
rocal of all nonzero Laplacian eigenvalues for U gi the next 
step is to find this sum, denoted by T g . By definition, 



N 



r * = E^ = E ^+ E (62) 

1=2 1 AfeAS, 1 ' 1 AfeA< 2) 1 



Let Tg 1 ^ and T g 2 ^ denote separately the two sums on the 



(63) 



rhsofEq. (J62j» . Obviously, 

4 1} = (/-!)(/ + I) 9 " 1 

(2) 

And T g can also be calculated as 



2(/+l) 




E 



1 t+U+i) 9 - 1 



i=2 \ A i+{ f + i)9-i A 1 + (/+1)B -i 



(64) 



Because A 9 and A 9 + ^ +1 j B _i are two roots of the quadratic 
equation given by Eq. (|3"21) . using Vieta's formulas, 
we have A 9 + A 9 +(/+1)s -! = Af -1 + / + 1 and A 9 x 
A 9 ,, .,,.<,_! = A 9-1 . Furthermore, considering A? = 0, 

i+\J + l) a 1 

so A 9 ,. , = m + 1. Then Eq. (l64l) is reduced to 
i+(/+i) a 

r (2) _ V- \ + / + 1 x + 



i=2 



= (/+ir i -i + (/ + i) 53 



i i 



A 9-1 / + 1 

i=2 J 



= (/ + l)f- 1 -l + (/ + l)T g _ 1 + 



/ + 1 



(65) 



Note that r g z> =T g - T g L> = T g -(f- 1)(/ + 
applying this result into Eq. (|62[) . one can reach the fol- 
lowing recursive relation for T g : 



r g = (f + i)r g _ 1 + f(f- 



1)8-1 f_ 



(66) 



With the initial situation To = 0, Eq. can be resolved 
to yield an explicit formula for T g as 



r g = (f + iy- 1 (fg-i) + 



Thus, the exact expression for F g is 



1 



/ + 1 



Fg = 2(f + iy- 1 (fg-l) + 



f + 1 



(67) 



(68) 
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which can be further represented as a function of network 
2/ 



size N g as 



F a = 



(/ + l)ln(/ + l) 



N g lnN g -^^N a + 2 



f + 1 y / + 1 



(69) 



When the network size tends to infinity, i.e., g — > oo, F g 
has the following dominant form 

F ^ (/+ i)tVi) Wg ' (70) 

a scaling also different from that previously found for 
Vicsek fractals^, in which F g increases as a superlinear 
function of N„. 



in the small-world trees, as shown above, it will visit at 
most a faction of nodes before being trapped; while 
for trapping in Vicsek fractals, the walker may visit a 
larger fraction (greater than j^j) of nodes prior to being 
absorbed by the central trap node. 



V. GENERALIZED GAUSSIAN STRUCTURES AND 
RELAXATION PATTERNS 

In this section, we consider the relaxation dynamics 
of the treelike polymer networks in the framework of 
GGS 3 »~— , which is an extension of the classic Rouse 
model^, developed for linear polymer chains and ex- 
tended to more complex geometries. 



C. Result comparison and analysis 



Brief introduction to GGS 



From above-obtained results given by Eqs. (|54[) 
and ([70|) . it is easy to see that the dominating terms 
for Fh (g) and F g behave differently. The former obeys 
F hg (g) ~ N g , while the latter follows F g ~ N g lnN g , 
greater than that of the former. This disparity indicates 
that in the family of treelike small-world polymer net- 
works, the location of the trap has a strong influence 
on the trapping efficiency measured by ATT, which is 
in comparison with that for Vicsek fractals, where the 
effect of trap's location is negligible 27 i 28 i 35 . In addition, 
the distinction between Fh (g) and F g also shows that 
the leading scaling of ATT to a given node in U g , e.g., 
the central node, might be not representative of the net- 
works. 

The dissimilar dominating scalings for Fh (g) and F g 
in U g lie in the network structure and can be heuristi- 
cally accounted for as follows. As shown in Fig. [2l U g 
consists of / + 1 copies of U g —i'. one central replica, and 
/ peripheral duplicates. When the trap is positioned at 
the central hub node, the particle will visit at most one 
copy of Ug-i, i.e., a faction of l/(/ + 1) among all nodes 
in U g . Thus, the ATT Fh (g) is small and grows linearly 
with network size, revealing a high trapping efficiency. 
In contrast, when the trap is located at another node, 
the particle should first visit the hub node, from which it 
continues to jump until being absorbed by the trap. So, 
the percentage of visited nodes is larger than that of the 
case when the trap is fixed at the hub. In particular, for 
the case that the trap is placed at a node farthest from 
the hub, the particle must visit all nodes of the networks 
before reaching the target. That is why the trapping 
process is less efficient when the trap is uniformly dis- 
tributed. 

The differences of behaviors of random walks in the 
small-world treelike polymer networks and Vicsek frac- 
tals are rooted in their underlying structures. For ex- 
ample, for trapping with a trap at the central node, the 
fact that the trapping efficiency of the former is higher 
than the latter can be understood as follows, for a walker 



The theory of GGS was accounted for in detail in pre- 
vious worka 3 ^— , thus we give here only a brief introduc- 
tion of the basic equation and main results related to the 
relaxation dynamics patterns. 

A GGS consists of N beads subject to the friction 
with friction constant £, which are connected to each 
other by elastic springs with elasticity constant K. In 
the Langevin formalism, the dynamics of bead m obey 
the following equation 



C^f^ + K Y, L ™i&m{t) = fm(t) + F m (t) . 



(71) 



In Eq. CO]), Rm(*) - (X m (t),Y m (t), Z m (t)) is the posi- 
tion vector of the mth bead at time t; L m i is the mith 
entry of the Laplacian matrix L describing the topol- 
ogy of the GGS; f m (i) is the thermal noise that is as- 
sumed to be Gaussian with zero mean value (f m (i)) and 
(fma(t)f m p{tf)) = 2k B T5 a p5{t - t'), where k B is the 
Boltzmann constant, T is the temperature, a and j3 rep- 
resent the x, y, and z directions; F m (t) is the external 
force acting on bead m. 

We focus on the motion (drift and stretching) of the 
GGS under a constant external force F = F@(t)e y (here 
0(t) is the Heaviside step function), switched on at t = 
and acting on a single bead in the y direction. The 
displacement along the y direction, Y(t), after averaging 
both over the fluctuating forces f m {t) and over all the 
beads in the GGS, reads^^ 



F 
aNC, 



N 1 

E 1 



exp(— aXit) 



A; 



(72) 



where a — K / '£ is the bond rate constant, and A, is the 
eigenvalues of matrix L with Ai being the unique least 
eigenvalue 0. 

Equation (|72p shows that, in the Rouse model the aver- 
age displacement depends on only the eigenvalues but not 
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the eigenvectors of matrix L. Notice that, in Eq. (1721) . 
due to Ai = 0, the motion of the center of mass has 
separated automatically from the rest. Moreover, from 
Eq. ([72)1 . the behavior of the averaged displacement for 
extremely short times and for very long times is obvi- 
ous. In the limit of very short times and sufficiently 
large N, (Y(t)) ~ Ft/(; while for very long times, we 
have (Y(t)) ~ Ft/(NQ. The physical explanation is as 
follows: for very short times only one bead is moving, 
whereas for very long times the whole GGS diffuses. The 
above two behaviors are general features for all systems, 
for a given GGS, its particular topology comes into play 
only in the intermediate time domain. 

In addition to (Y(t)), another interesting quantity is 
the mechanical relaxation form, namely the complex dy- 
namic modulus G * (uj), or equivalently, its real G'(w) 
and imaginary G"(uj) components, which are known as 
the storage and the loss modul i 41 ' 42 . For very dilute so- 
lutions and for uj > 0, G'{ui) and G"(oj) for the Rouse 
model are given by 



N 



N 



and 



^ 1 + (uj/2a\ t ) 2 



N 



„ vk B T ^ co/2aX l 

N ^ 1 + ( w /2aA,) 2 



(73) 



(74) 



where v denotes the number of polymer segments (beads) 
per unit volume. 

The relaxation patterns of various polymer systems 
have been studied in previous works 1 -, including star 
polymer s 38 ' 39 , dendrimer o 14 ' 43 " — , hyperbranched poly- 
mers 4 ^—, dual Sierpinski fractals 5 ^— , small-world net- 
works 54,55 , and scale-free networks 5 ^. Below will com- 
pute related relaxation quantities for the treelike small- 
world polymer networks under consideration. 



B. Relaxation patterns 

By substituting the full eigenvalues obtained in sec- 
tion |nrj] into Eqs. ([72"]). (1731 . and ([711), we can compute, 
respectively, the averaged displacement (Y(t)), the stor- 
age modulus G'(oj) and the loss modulus G"(u>) for the 
relaxation dynamics of the small- world polymer networks 
U g . 

We begin by focusing on the averaged monomer dis- 
placement, (Y(t)), given by Eq. ([72")) in which we set 
a = 1 and F/( = 1. In Fig. [5] we present in a dou- 
ble logarithmical scale the results of (Y(t)) for networks 
Uq with / ranging from 2 to 6. As mentioned above, 
from Fig. [5l the behavior of (Y(t)) for very short and 
long times are clearly evident, obeying (Y(i)) ~ Ft/C, 
and (Y(t)) ~ Ft/(NQ, respectively. In the region of 
very short times, only one monomer moves, hence the 
curves are not dependent on N. In contrast, in the 



domain of very large times, the whole structure drifts, 
thus the curves depend on N: the higher the value of 
N, the slower the limiting long time behavior will be. 
Typical for the small- world treelike structure is interme- 
diate time regime, where (Y(t)) scales as a power-law 
behavior with the exponent a = 0.2 for all /, a phe- 
nomenon different from that of Vicsek fractals, the ex- 
ponent of which is related to their spectral dimensions 
d=21n(/+l)/ln(3/ + 3). 




icr 2 io° io 2 io 4 io 6 io s io'° io 12 
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FIG. 5. (Color online) Averaged monomer displacement 
(Y(t)) for U 6 with / = 2,3,4,6. 

For the storage modulus G'(uj), we report the results in 
Fig. O which is plotted in dimensionless units by setting 
(7 = 1 and vk ^ T — 1. Figure [6] indicates that in the 



N 

very low and high frequency limit the storage modulus 
G'{uj) exhibit a power-law uj 2 and a plateau, respectively. 
Both phenomena are the same as those of many different 
systems. In the intermediate regime the structure being 
studied play an important role. For the four cases of / = 
2,3,4,6, we can observe an obvious power-law behavior 
with an exponent a' = 1 for all /, but the behavior 
becomes more prominent with / increasing from 2 to 6. 
It is worth stressing that this result is also different from 
that for Vicsek fractals^—. 

For the loss modulus G"(u>), we plot in a double scale 
the results in Fig. [7J As in the case of G'(u>), we consider 
fj = l and §r T — 1. From Fig. [TJ it is easy to notice 
that for very low frequencies w, G"(uS) ~ uj 1 ; and that for 
very high frequencies ui, G"(ui) behaves as G"(ui) ~ u^ 1 . 
In the intermediate region, no power-law behavior is ob- 
served, which is in marked contrast to that corresponding 
to Vicsek fractals 4 ^— . It is also important to notice that 
in the intermediate region, G'(ui) and G"{lS) display dif- 
ferent behavior for the small-world structure. 

The distinct behaviors for the three quaternities re- 
lated to relaxation patterns in Viscek fractals and the 
small-world treelike polymer networks lie in the differ- 
ences between the two structures. As the name suggests, 
Viscek fractals are fractals, their relaxation patterns are 
determine by the fractal dimension and spectral dimen- 
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FIG. 6. (Color online) Storage modulus G'(oj) for Ue with 
various /. 




FIG. 7. (Color online) Loss modulus G"(u>) for Ue with / = 



sion^Zr— . For the small-world treelike polymer networks, 
they are non-fractal, and thus exhibit different relaxation 
patterns. 



VI. FLUORESCENCE DEPOLARIZATION 

We are now in position to study the dynamics 
of Forster energy transfer over a system of chro- 
mophore o 2 ' 3 ' 15 positioned at nodes (beads) of the small- 
world polymer networks. We suppose that the energy can 
be exchanged only between the nearest neighbors. Then, 
the energy transfer among chromophores located at the 
nodes of U g can be described by the following equation 



dPi{t) 



A'., 



Pi(t), (75) 



where Pi (t) denotes the probability that node i is excited 
at time t and represents the transfer rate from node 
j to node i. 

As usual, we here separate the radiative delay (equal 
for all chromophores) from the transfer problem. In fact, 
the radiative delay only leads to the multiplication of 
all the Pi(t) by exp(— g/ru), where 1/tr is the radiative 
decay rate. We presume that all microscopic rates are 
equal to each other, say k, then Eq. [75] becomes 



dPjjt) 
dt 



Kg 

-h^ujp 



'j(t) - (kLi^j Pi(t) , (76) 



where Lfj is the ijth entry of Laplacian matrix L g . 

As shown before^i 3 -^, the probability of finding the ex- 
citation at time t on the originally excited chromophore, 
averaged over all possible starting points on U g , is given 
by 



which is dependent on all eigenvalues of the Laplacian 
matrix for U g . 

Making use of the eigenvalues obtained in Section lLU Al 
we can evaluate (P(t)) for very large networks, without 
diagonalizing the Laplacian matrix. By setting k = 1, 
i.e., by measuring the time in units of 1/fc, we can com- 
pute the average probability (P(t)) that an initially ex- 
cited chromophore is excited at time t. In Fig. [51 we 
present the results for the case / — 3, with g varying 
from g — 4 to g = 7. 

From Fig. [8j we can see that at very short and very 
long times, the overall behavior for different g is similar. 
For example, at long times (depending on the network 
size), each curve becomes flat, which (in the absence of 
any radiative decay) is due to the equal distribution of 
the energy over all nodes in the networks, with each node 
having a probability of l/N g of being excited. We note 
that similar phenomenon is also observed for Vicsek frac- 
tals^ 3 -. However, at intermediate times, the curves for 
different g behave quite different, but no scaling is ob- 
served, meaning that no curves follow a linear behavior. 
This phenomenon is as opposed to that for Vicsek fratals, 
the corresponding curves of which show an obvious alge- 
braic behavior^ 3 -. The disparity in (P(t)) makes it easy 
to differentiate between Vicsek fractals and the polymer 
networks studied here. 



VII. CONCLUSIONS 

In this paper, we have introduced a class of dctcrmin- 
istically growing treelike polymer networks, and shown 
that they have an exponential-form degree distribution 
and the small-world characteristic at the same time. We 
have fully characterized the Laplacian eigenvalues and 
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FIG. 8. (Color online) The average probability (P(t)) for 
/ = 3 and g =4, 5, 6, and 7 from above, shown in a log-log 
scale. 



their corresponding eigenvectors of the networks, which 
are determined through recursive relations derived from 
the specific network construction. Using the eigenvalues, 
we have further studied three representative dynamics for 
the polymer networks, such as trapping problem, relax- 
ation dynamics in the framework of the GSS, and energy 
transfer through fluorescence depolarization. Moreover, 
we have compared the dynamical behaviors with those for 
Vicsek fractals, which are fundamentally different from 
each other. Finally, in addition to the aforementioned 
dynamics, we expect that the obtained eigenvalues and 
eigenvectors can be adaptable to other dynamics in the 
small- world networks, e.g., quantum walks^Zr— . 
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